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We propose a generalization of the Bjorken in-out Ansatz for fluid 
trajectories which, when applied to the (1+1) hydrodynamic equa- 
tions, generates a one-parameter family of analytic solutions interpo- 
lating between the boost-invariant Bjorken picture and the non boost- 
invariant one by Landau. This parameter characterises the proper- 
time scale when the fluid velocities approach the in-out Ansatz. We 
discuss the resulting rapidity distribution of entropy for various freeze- 
out conditions and compare it with the original Bjorken and Landau 
results. 

1 Introduction 

There is an accumulating evidence that hydrodynamics may be relevant for 
the description of the medium created in high-energy heavy ion collisions 
PQ. Indeed, experimental measurements such as the elliptic flow [2] shows 
the existence of a collective effect on the produced particles which can be 
described in terms of a motion of the fluid. More precisely, numerical sim- 
ulations of the hydrodynamic equations describe quite well the distribution 
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of low-p± particles [T], with an equation of state close to that of a "perfect 
fluid" with a rather low viscosity. This evidence is of course indirect, since 
it relies on assumptions about the initial and final stages of the evolution of 
the fluid. Thus some doubts can be cast either on the full thermalization 
of the medium, or on the possibility of accounting for some viscosity of the 
fluid [3 J . Also, hydrodynamics are not expected to work for leading particles 
i.e. near the kinematic light-cone. 

Given these objections, it is important to separate precisely the conse- 
quences of the hydrodynamic flow from those of the initial and final con- 
ditions. From that point of view, it seems tempting to discuss a simplified 
picture which can be qualitatively understood in physical terms. One such 
simplification, which we are going to follow in this paper, is the idea that 
the evolution of the system before freeze-out is dominated by the longitudi- 
nal motion [U [5] and thus, in fact, the hydrodynamic transverse motion can 
be neglected or at least factorized out. Thus we shall consider the (1 + 1) 
dimensional system. 

On the theoretical grounds, there are also quite appealing features for 
applying hydrodynamic concepts to high-energy heavy-ion reactions. Such 
concepts have been already introduced some time ago and find a plausible 
realization nowadays. The fact that a rather dense medium is created in the 
first stage of the collision allows one to admit that the individual partonic or 
hadronic degrees of freedom are not relevant during the early evolution of the 
medium and justifies its treatment as a fluid. Moreover, the high quantum 
occupation numbers allow one to use a classical picture and to assume that 
the "pieces of fluid" may follow quasi-classical trajectories in space-time, 
expressed as an in-out cascade |6j with the straight-line trajectories starting 
at the origin, i.e. 

y = v (i) 

where 

1 , / E + p\ 1 (t + z\ 

are respectively the rapidity and "space-time rapidity" of the piece of the 
fluid. 

Note, for further use, that ([1]) can be rewritten in the form 

2y = logw + — logw~ = \ogz + — log;z~ (3) 

where w ± = e ±y are the light-cone components of the fluid (four-)velocity 
and z± = t ± z are the light-cone kinematical variables. 
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Taking (pQ) as the starting point and using the perfect fluid hydrodynam- 
ics, Bjorken developped in his seminal paper [5] a suggestive (and very useful 
in many applications) physical picture of the central rapidity region of highly 
relativistic collisions of heavy ions. In this picture the condition ([1]) leads to 
a boost-invariant geometry of the expanding fluid and thus to the central 
plateau in the distribution of particles. 

It is now experimentally established [7] , however, that the central rapidity 
region of heavy ion collisions is better described by a Gaussian fit with a width 
proportional to Y = logs, the total rapidity range of the secondaries. This 
finding has renewing interest [8] for the pioneering hydrodynamic description 
by Landau [I] where, indeed, a Gaussian-like distribution of the fluid was 
obtained. For the same reason, new families of 1+1 relativistic hydrodynamic 
solutions have been recently proposed 0, [10] ■ 

In the present paper we propose to study a generalization of the formula 
(JTJ) for the classical trajectory which (as we show in the following) interpolates 
naturally between the Landau and Bjorken pictures: 

2y = \ogu + - logu" = log /+(«+) - log /_(«_) (4) 

where f±(z±) are a priori arbitrary functions. They have to be determined 
from the hydrodynamic equations. 

The hydrodynamic equations are rewritten in terms of light- cone vari- 
ables in the next section. The consequences of the quasi-classical Ansatz 
(JTJ) and of the generalized one (jlj) are dicussed in Section 3 where also the 
corresponding solutions of the hydrodynamical equations are derived. Var- 
ious selections of the freeze-out conditions are discussed in Section 4. Our 
conclusions and comments are listed in the last section. 



2 Hydrodynamic equations in light-cone vari- 
ables 

We consider the perfect fluid for which the energy-momentum tensor is 

= (e + p)u"u v -prf v (5) 

where e is the energy density, p is the pressure and is the 4-velocity. We 
assume that the energy density and pressure are related by the equation of 
state: 

e = gp (6) 
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where 1 / y/g is the sound velocity in the liquid. Using 

u ± = u°±u 1 = e ±y (7) 

and introducing 

2± = i±2 = ,« ± -=^ (^±^)= ^ Jfe (8) 

where r = ^/z + z_ is the proper time and rj is the spatial rapidity of the fluid, 
the hydrodynamic equations 

dpT*" = (9) 

take the form 

d±T m + -d+(T n ± T 00 ) - -a_(T n =f T 00 ) = . (10) 
Using now ((SJ) and the equation of state ([6]) we deduce from this 

#d+[logp] = — d+y — e 

00-pogp] = ^l!a_y + ^e 2 ^ + y . (11) 

These are two equations for two unknowns which describe the state of the 
liquid: the pressure p and the rapidity y. They should be expressed in terms 
of the positions z + , z_ in the liquid. Other thermodynamic quantities can be 
obtained from the equation of state (jSj) and the standard thermodynamical 
identities: 

p + e = Ts; de = Tds (12) 

where we have assumed for simplicity vanishing chemical potential. 
The result is 

e = gp = eo T 9+1 ; s = s T 9 -> s ~ e a/{9+1) . (13) 
Note that (TTTT) implies the consistency condition 

d + d-y = -^y 2 {d^[e- 2y } - d + d + [e + **\} . (14) 
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3 Generalized in-out Ansatz 



3.1 Bjorken's in-out Ansatz and boost- invar iance 

The simplest possibility to describe the expansion of the fluid was suggested 
by Bjorken |5] who proposed to use the Ansatz ([T]) in the hydrodynamical 
context. Introducing into ffTTj) we obtain 

gd + [\ogp] = ! 00_[logp] = ( 15 ) 

from which we deduce 

p = e g- 1 = p (z + z_)- {9+1 ^ 29 = p T ~^ /9 , (16) 

where po is a constant. Thus the system is boost-invariant: the pressure does 
not depend neither on rj nor on y. So are e, s and T, given by ( 1131) . 

3.2 Beyond boost invar iance 

The data on both nucleon-nucleon and nucleus- nucleus collisions (see, e.g. 
[7]) show that the produced system strongly violates boost invariance (except 
perhaps in a narrow region of small cm. rapidities). It is thus necessary to 
go beyond ([I]). As already indicated in the Introduction, we propose to study 
-as a simplest generalization of ([T|)- the Ansatz (jlj). Introducing (jl]) into (fT4j) 
we obtain 

f-d_d_(f-) = Ud + d + (f + ) = A 2 /2 (17) 
where A 2 is a constant. Thus both /+ and /_ satisfy an identical equation: 

ff" = A 2 ! /2. (18) 

Note that A = implies /" = and thus we recover the condition (JTJ). 
We conclude that A describes the deviation of the system from the uniform 
Hubble-like expansion. 

Eq. (ITS1) can be solved multiplying by /' and dividing by /: 

[(f) 2 ]' = A 2 [\ogf]' - /' = Ay/log(f/H) (19) 
where H is an arbitrary constant. 
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Eq. (1191) can be solved in the standard manner. We obtain 

dF' 

7Hp 



, F dF' 

z- z = h I — (20) 

F 



where we have introduced the notation 

F = f/H- h — H/A . (21) 
When (j4|) is introduced into f|T9|) we obtain 

{l+gfF' + g*-lF>_ 



gd + \logp\ 



4 F, 



gd„[\ogp] = -LJ^^ + *L__±. ( 22 ) 
From this we deduce 

glogp = — logF + + ^—F_ j — + A_(z_) 

glogp = -^±^logF_ + J ^ + A + (z_) . (23) 

The integrals on the R.H.S. can be evaluated using ( |T9l) . Indeed 

/•<fz 1 , f d\ogF n1 A . . 

— = — — = h / - 6 = 2/iv/logF . 24 

J F J F F> J v4oiF V 6 1 ; 



We thus obtain 



glogp = log(F + F_) + ^v/log(F + )log(F_ 



(25) 



where, for the two equations (1231) to be consistent with each other, we had 
to take A±(z±) = — lH ? logFj-. This finally gives 

p(z + ,z.)=p expj-^^^ + q+^M-} (26) 

V(*H *-) = ^ 0+ - il) (27) 

6 




0.4 
0.2 

I , , , , , 

20 40 60 80 100 120 

x/h 

Figure 1: The ratio y/rj plotted versus proper-time r/h for several values of 
space-time rapidity tj. The boost-invariant picture corresponds to y/rj = 1. 
The "full-stopping" initial condition is y/rj = at r = 0. The asymptotic- 
value is still 1, but it is reached only at very large proper-times. 



where, by definition 

l±(z±) = y/log(F±) . (28) 

This completes the solution of hydrodynamic equations constrained by 
the generalized in-out cascade (jlj). One sees that the pressure depends on 
both r and rj. Thus the system is not boost-invariant but boost invariance is 
recovered in the limit h — > 0, z± fixed (see Eq. (1201) ). On the other hand, it 
can be remarked that the Landau asymptotic solution [1] can be recovered 
in the limit h fixed, z± — *■ oo. 

Other thermodynamic parameters are obtained from ffl3|) . giving 

8 ~ <*/<**) = s exp {-^p K + l -} + ^ M-} = exp (- 9 0) (29) 
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where we have denoted 

9 = log(T /T) = ±±l[ll + ll] - izl i+L (30) 

with s and T denoting the entropy and temperature at the beginning of 
the hydrodynamic evolution. 

To illustrate the deviation of our solution (|27p from the in-out Bjorken- 
Gottfried-Low Ansatz (jTj), we show in the Figfj] the ratio y/rj vs r/h (the 
proper time measured in units of h) for several (fixed) values of 77. One sees 
that in the limit r/h — > the rapidity y vanishes. In this limit we have (e.g. 
for the region y > 0) 

y « (r/2/i) 2 sinh 2r/ = (t/2h) 2 - (z/2h) 2 < ^(t/hf . (31) 

Thus for a fixed (small) t/h the fluid starts at rest and acquires some 
velocity at later times, as in the Landau "full stopping" solution. At large 
times, t — > oo, one obtains y ~ i], i.e. the in-out Ansatz ([T]) is approximately 
recovered. Thus our solution does indeed interpolate between the Landau 
and Bjorken hydrodynamics. 

A last comment is in order. In all cases, the solution of the flow is also 
defined outside the kinematical light-cone. Indeed, there is some flow of 
energy entering the light-cone from outside. It could be physically inter- 
preted as mimicking energy sources on the light-cone ("leading particle ef- 
fect"). However, the relevance of hydrodynamical models near the light-cone 
is questionable. 

4 Entropy density at freeze-out 

The observable results of the model depend in an essential way on the as- 
sumed freeze-out surface. The point is that the densities s and e which enter 
the hydrodynamic equations are densities per unit volume in the rest frame 
of the fluid. But we are generally interested in the distribution of entropy 
dS/dy and/or of energy dE/dy densities per unit of rapidity, as these quan- 
tities are possible to measure. For given s and e, dS/dy and dE/dy depend 
on the hypersurface at which the hydrodynamic evolution stops and the fluid 
changes into particles (freeze-out surface). To fix attention, in the following 
we discuss the entropy density. 



8 



4.1 General freeze-out surface 

The evaluation of the entropy density per unit of rapidity for a given freeze- 
out surface can be performed in two steps. 

First we evaluate the amount of entropy in an infinitesimal volume along 
the freeze-out surface: 

dS = su^da^ (32) 

where is the 4-velocity of the fluid and da^ is the 4-vector orthogonal to 
the surface satisfying 

da^da^ = dz^dz^ = dz + dz_ (33) 
Consider the (space-like) surface 

z_)=C (34) 

where C is a constant. Then 

<$> + dz + + = ; $ ± = — — . (35) 

oz± 

It follows that the unit vector orthogonal to the surface is 

n+ = 7d= ; n - = T^W- ■ (36) 

The infinitesimal length along the surface is 

da = ^Jdz + dz_ = dz_ a/$_/$ + = dz + a/$ + /$_ (37) 

Therefore 

u»da» = [u^ + + u-^}^- = [u + * + + rt-]^ • (38) 

In the second step we express the infinitesimal volume along the freeze- 
out surface in terms of the infinitesimal interval of rapidity. This can be done 
using the relation (TSUj) which gives z± as a function of F± = exp (/|) . We 
have 

dz± dF± 2 exp (lp 
dZ± = dF±M[ = k dl ± (39) 
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Using this relation and (1201 . the R.H.S. of (1551) can be expessed in terms of 
l + and /_. 

This in turn can be expressed in terms of rapidity y using (j3J) and the 
condition (134|) which gives an additional relation between z + and z_ and thus 
following (120]) also between l + and /_. In particular, using the differential 
forms, we have 

„ 2 _ ±2$ T e^// T _ ±2fe$ T 

^ $ + exp(^)// + +$_exp(^)//_^ "> d2± " « &+ L e -£ + *_* +e V y 

(40) 

and thus, finally 

4.2 Freeze-out at a fixed time 

In this section we take the surface at t — const, for a first example. In the 
notation from the previous section we write 

=t= ^(z + + z_) = const; $+ = $_ = - (42) 



Using (T4T|) we thus have 



Ml+i 2 -) 



s ea v + 

dfS 1 = su^do u = — — dy = 

M /+ + /_ + (/ + -/_)tanhy y 

Sq e -(^l)(i+-i-) 2 /4 

/+ + /_ + (Z+ - Z_)tanhy ^ ' ^ 

If, following Landau, we approximate both l + and /_ by large constants, 
then for finite y the difference (l + — Z_) is small and we have 

~ e -(^)('+-'-) 2 /4 . (44) 

For (7 = 3 this formula is identical to the asymptotic result of Landau 
pEj. This can be seen when displaying the distribution dS/dy; In FigfS], one 
shows dS/dy with t = const, (formula (|43l ) compared with the Landau 
approximation, formula (jSJ), for different values of the parameter h, which, 
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Figure 2: The curve dS/dy with t = const, (cf. formula (|43|) ) compared 
with the Landau approximation (cf. flSJ), dashed lines). The kinematical 
end-points at y ma x correspond to 2_ = (y > 0). 

by simple rescaling of the kinematic variables, correspond to different end- 
points in rapidity. 

One should keep in mind, however, that the relation t = const, implies a 
different relation between l± and y than the condition r ~ const., which is 
the freeze-out condition considered.] by Landau [I]. As discussed in the next 
subsection, this leads to a rather different shape of the distribution dS/dy. 

4.3 Freeze-out at a fixed proper time 

To investigate the relation to the Landau solution and its comparison with 
the Bjorken one we consider the freeze-out at a fixed proper time. 

^^More precisely, Landau discusses the limitation of the 1+1 dimensional motion by its 
transition to the 3+1 dimensional hydrodynamical motion and shows that it appears at 
t ~ const, playing the role of a freeze-out surface. 
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Using the notation of Section 5.1 we have 

§(z + ,z_) = z + z_ = t 2 = const. (45) 
z + dz^ + z_dz + = ; $ ± = z T (46) 



and thus 



dS = he-^ l +- 1 ^ 2 " eVz - + e ~ Vz + dy . (47) 

LeVz_ + l+e~yz+ 



This is a general formula. When supplemented by (j4]) and (120]) . it expresses 
dS in terms y and r. 

When /i-^0we can use the approximation (see the Appendix) 

z± = (48) 

to obtain 

rf5 = /, !± + lz e -^a + -^-) 2 rfy . (49) 
.£t_t_t_ 



For Z-|_ — >■ oo and fixed y one recovers the Landau result 

The result given by (fi"7j) is plotted in FigJ3] where dS/dy, is displayed 
for different values of the parameter h and compared with the approximate 
formula 

^ = So ev 7 ^ (50) 
dy 

obtained by Landau jl]. The parameter L was adjusted to obtain the correct 
slope at y — 0. 

One observes some deviations from the perfect Gaussian which was con- 
sidered in a simplified version [12] of the Landau model (and which agrees 
-if the multiplicity distribution is assumed proportional to the entropy- with 
the data [7]). Note that at fixed r and h — > 0, the distribution becomes 
significantly flatter, going smoothly to the Bjorken limit at h = 0. 



4.4 Freeze-out at a fixed temperature 

Instead of considering the freeze-out surface at the limit where the transverse 
motion becomes relevant (cf. [1]), a natural conjecture is to consider freeze- 
out at a fixed temperature, i.e. when the temperature reaches the value 
where pions are expected to become liberated, e.g. [H ITT]. 
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Figure 3: The curve dS/dy with r = const, for various values of r/h. Full 
lines: Eq. ( I49|) ; Dashed lines: Eq. (|50p . Note that, for clarity, the curves were 
shifted by log 10 r/h. 



Denoting the initial temperature by T and the freeze-out temperature 
by Tp we have 

z.) = ^ \ll + l + l_ = g log(T /T F ) = g e = const. (51) 



Hence 



= ^[(g+l)l+-(g-l)l-]l+ = 2 



2 LU? 1 y + v ^ > ~' + 4^exp(/2 
where we have used the relation following from ( 120|) : 
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Therefore using (1251) and (14T1) we have 



( 5 + l) M __£i( Z 2 +Z 2^ • 



Now, the relations 
1 + 5 



R + q 



5-1 



l+l-=g9;; ll-£ = 2y 



imply 



giving 



(54) 

(55) 
(56) 



(g + !)/+/_ 



2 

5-1 



2g^9 2 -y 2 /g; 



9-1 



- ^Q2 _ y2/g 



i + + i_ = sHy [e - v^ 2 - y 2 /g 

Thus we finally obtain 

dS ~ e^k^M (e - v^TT^) 



-1/2 



-1/2 



ydy 



(57) 



(58) 



v^ 2 - y 2 /g ' 

One sees that this formula exibits a singularity by its transition to the 
3+1 dimensional hydrodynamical motion at y 2 = g6 2 which is of course 
unphysical and reflects the singular initial conditions of our approach and the 
expected limitation of hydrodynamic models to the more central scattering 
region. 

It is, however, interesting to observe that the hypersurface T =const is 
only partly space-like. It becomes time-like at the rapidity determined from 
the condition $_ = 0, i.e., (g + 1)Z_ = (g — l)l + , giving (c.f. (|55|) ) 



Un 



2g 
5+1 



0; l 2 + = {g + l)9. 



(59) 



This is illustrated in FigJH where two profiles 9 = const, are displayed. 



2 The singularity in dS/dy does not come from a true singularity in the kincmatical 
domain but is due to a maximum value reached by the rapidity y as a function of z_. 



14 



FIXED TEMPERATURE PROFILES 

x = const 




1 2 3 4 5 6 z 

Figure 4: The profiles 9 = const, for two values of 9 = log T /T (continuous 
lines). The comparison is made with fixed r (dashed lines). Note that a 
different value of the parameter h has been chosen to obtain a comparable 
range in space-time. 

The numerical estimates show that -for large enough 9 (9 > 2)- the effect 
of the singularity on the entropy distribution as a function of y is not signifi- 
cant in the region y < y m . This is shown in Figj5] where dS/dy is plotted for 
several values of 9 > 2. One sees that all distributions are close to Gaussians. 
The slope, however, is rather small, certainly smaller than required by data, 
unless one considers a larger value of the parameter g, i.e. a smaller value of 
the speed of sound. 



5 Conclusions and outlook 

We have investigated longitudinal hydrodynamic expansion of a perfect fluid 
forming an infinitely thin layer at the initial time and satisfying the equation 
of state with an arbitrary sound velocity. We proposed a generalized in- 
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Figure 5: The curve dS/dy with 8 = const, for various values of 8, y < 38/2. 



out Ansatz which unifies the Landau and Bjorken approaches to relativistic 
hydrodynamics. 

Our conclusions can be summarized as follows. 

(i) When the Bjorken-Gottfried-Low in-out hypothesis (JTJ is assumed, 
the hydrodynamics implies that all thermodynamic properties of the fluid 
depend only on proper time, i.e. the system is boost-invariant. 

(ii) We proposed a generalized in-out hypothesis, (@J, and discussed its 
implications. It turns out that in this case the hydrodynamic equations can 
be solved in an analytic form, giving explicit formulae for the thermody- 
namic characteristics of the fluid in terms of their initial values and one free 
parameter h defining a dynamical scale in configuration space. 

(iii) The resulting entropy distribution in rapidity, dS/dy, was evaluated 
and shown to depend significantly on the assumed condition for the freeze- 
out. 

(a) For freeze-out at a fixed proper-time the density is close to (but with 
some deviation, particularly at small r, from) a Gaussian which is tradi- 
tionally attributed to Landau solution. It tends smoothly to the Bjorken 
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boost-invariant solution for h — > 0. 

(b) For freeze-out at a fixed temperature the distribution is even closer 
to a Gaussian, if one restricts to the region where the freeze-out surface is 
space like and if the ratio T /Tp is large enough. 

(iv) It is worthwhile to note that the freedom in the choice of the value of 
the sound velicity may be helping in phenomenological applications of these 
results to data. 

Compared to other recent (1 + 1) hydrodynamical models [9j [ID] , our 
solution is mainly characterized by the smooth, one-parameter dependent 
transition between the Bjorken and Landau hydrodynamical models and by 
its analytic simplicity. It would be useful to study further the classification 
of all the solutions in a unified framework. 

More generally, there is a clear need for an extension of our investigation 
to include more flexible initial conditions, relaxing the point-like character of 
the fluid at the beginning of the evolution. This, however, demands a more 
sophisticated analysis (e.g. an application of the general recipe of [13]) and 
goes beyond the scope of this paper. 

On a theoretical ground, it would be interesting to have a physical in- 
terpretation of the generalized Ansatz (j3J), which appears as a mathematical 
harmonic property d + d-y = of the hydrodynamical flow. In particular, 
an extension to this case of the Gauge/Gravity correspondence applied in 
Ref. [14J to the Bjorken Boost-invariant flow, would be insightful. 
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Appendix. Solution of the equation (1201) 

We rewrite (1201) as 



{z-Q = h 



dv 



F 



(60) 



where F = f/H and h = H/A. 

Changing the variable of integration: 



logv = u ; 2udu = dv/v ; dv/u = 2e u du 



(61) 



we arrive at 



z - C = 2hF 



D (y^gFj - D (VlogFo) 



(62) 

This integral is the so called Dawson's integral: -D(x) = e x J e M du. For 
large A/log F it approaches (2yTog F) _1 and thus we obtain 



Av^ogF 



(63) 



The asymptotic expansion of D(x) is 



1 ; 2s A? r(l/2) x 2 " 



n=0 



(64) 



For small a; one can evaluate this integral effectively by the series expan- 
sion: 



D(x) = e~ x 
Thus we have in this case 



E 

n=0 



x n 



00 2n 
,2 X M 



xe 



nl 



E 

n=0 



(2n+l)n! 



(65) 



H 

z - z = 2— 



A TEn (log^T /T^TV- (logFo)" 

n=0 v 7 n=0 v 1 



(66) 
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